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We propose a practical scheme for calculating the local gravitational self-force experienced by 
a test mass particle moving in a black hole spacetime. The method — equally effective for either 
weak or strong field orbits — employs the mode-sum regularization scheme previously developed for 
a scalar toy model. The starting point for the calculation, in this approach, is the formal expression 
for the regularized self-force derived by Mino et al. (and, independently, by Quinn and Wald), 
which involves a worldline integral over the tail part of the retarded Green's function. This force is 
decomposed into multipole (tensor harmonic) modes, whose sum is subjected to a carefully designed 
regularization procedure. This procedure involves an analytic derivation of certain "regularization 
parameters" by means of a local analysis of the Green's function. This manuscript contains the 
following main parts: (1) Introduction of the mode sum scheme as applied to the gravitational case. 
(2) Two simple cases studied: the test case of a static particle in flat spacetime, and the case of a 
particle at a turning point of a radial geodesic in Schwarzschild spacetime. In both cases we derive 
all necessary regularization parameters. (3) An Analytic foundation is set for applying the scheme in 
more general cases. (In this paper, the mode sum scheme is formulated within the harmonic gauge. 
The implementation of the scheme in other gauges shall be discussed in a separate, forthcoming 
paper.) 
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I. INTRODUCTION 

A pointlike object of mass m (hereafter a "particle") is known to move along a geodesic of the background spacetime 
in the limit where m — s- 0. If the mass of the particle is finite, the motion would no longer be geodesic: interaction of 
the particle with its own gravitational field will then give rise to the exertion of a "self force" , pushing the particle away 
from geodesic motion. (The non-conservative part of this force is customarily referred to as "radiation reaction".) 
The problem of calculating the gravitational self force is a longstanding one. This problem is usually tackled in 
the context of perturbation theory, by treating m as a small parameter and looking for the 0{m) correction to the 
geodesic equation of motion on a fixed background. The prototype problem — calculating the electromagnetic self force 
on an accelerating (classical) electron in flat space — was considered by Dirac in his famous 1938 paper [Q. Already in 
this flat space problem one encounters the fundamental issue of regularization, namely, how to correctly handle the 
divergence of the electromagnetic field (and self force) at the particle's location. Dirac's regularization yielded what 
is now referred to as the "Abraham-Lorentz-Dirac" (ALD) force (proportional to the time derivative of the electron's 
acceleration). The interpretation of the ALD equation of motion and its solutions has been subject to much further 
study 1^ . A formal framework for calculation of the electromagnetic self force in curved spacetime was first developed 
by DeWitt and Brehme in 1960 §. Here, in addition to re-addressing the question of regularization, one must also 
deal with the non-local nature of the self force effect: In curved spacetime, waves emitted by the particle at a given 
instant may backscatter off spacetime curvature and interact back with the particle later on its motion. Consequently, 
the momentary self force acting on the particle appears to depend on the particle's entire past history. 

The continuous effort for theoretical characterization of gravitational waveforms from strongly gravitating astro- 
physical systems led, in recent years, to a rise of interest in the problem of calculating self forces in curved spacetime. 
The most prominent consequence of this has been the first development of formal tools for calculating the gravita- 
tional self force (see below). Knowing the self force is necessary for describing the strong field orbital evolution in 
various astrophysical scenarios, in particular — the capture of a small compact object by a supermassive black hole 
[an event expected to serve as a main source of gravitational waves for the future space-based Laser Interferometer 
Space Antenna (LISA) In some occasions, the (time-averaged) orbital evolution, and the consequent gravitational 
waveforms, may be determined by calculating energy-momentum fluxes at infinity (and through the event horizon) 
and using balance considerations [||. This method can be applied successfully in models where the central mas- 
sive object is spherically symmetric (e.g., a Schwarzschild black hole), or even in the more realistic case of a Kerr 
background — but then, only for equatorial orbits. In more general cases (a non-equatorial orbit around a Kerr black 
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hole), a full specification of the orbital evolution requires knowledge of the temporal rate of change of the Carter 
"constant of motion" , which, to the best of our knowledge, cannot be achieved by standard balance considerations. 
In addition, balance considerations involve averaging over orbital periods, and can only account for situations where 
the orbital evolution is adiabatic. This method is inadequate in other, non-adiabatic scenarios, as the final plunge of 
the particle into a black hole. In general, therefore, it seems that one cannot avoid tackling the problem of calculating 
the local, momentary self force. 

In 1997, two groups— Mino Sasaki and Tanaka (MST) ^ and Quinn and Wald (QW) 0— have worked out (in- 
dependently) a formal framework for local gravitational self force calculations in curved spacetime. Their combined 
work presents three different derivations of the self force, all yielding the same formal expression (though in a different 
level of rigorousness). In the first approach, MST ||] directly generalized the above mentioned analysis by DeWitt and 
Brehme to the gravitational case: First, the gravitational perturbation induced by the mass particle was evaluated 
near the worldlinc using Hadamard's expansion of the retarded Green's function Then, the equation of motion was 
deduced by imposing local energy-momentum conservation on a thin worldtube surrounding the particle's worldline. 
The second approach, still by MST employed the technique of matched asymptotic expansion, which is based on 
considering two asymptotic zones outside the particle: an "internal" zone, where the geometry is taken to be that of a 
Schwarzschild black hole with a tidal perturbation associated with the background curvature, and an "external" zone, 
where geometry is that of a perturbation on a fixed background. The equation of motion is then obtained by requiring 
consistency of the metric in the matching region of the two zones. The third approach for deriving the gravitational 
self force was presented by QW (7|, and is based on what they called the comparison axiom. According to this axiom, 
the correct self force can be deduced by appropriately comparing the perturbation on the given curved spacetime with 
that on flat spacetime. This procedure results in the elimination of the divergent piece of the force, and, presumably, 
in extraction of its correct finite part. (QW also applied their approach for calculating the electromagnetic |7j and 
scalar Q self forces.) 

As already mentioned, all three derivations of the gravitational self force yield the same result, which we now briefly 
describe. In curved spacetime, the two-point retarded Green's function associated with a certain wave equation (i.e., 
the gravitational perturbation equation in a given gauge) is composed of two parts: an "instantaneous" part, which 
describes the propagation of influence along the light cone; and the so-called "tail" part, describing the nonlocal effect 
of waves propagating inside the light cone. The gravitational self force is derived from the gravitational perturbation 
produced by the particle, which, it turn, can be expressed as a worldline integral over the Green's function. MST 
and QW (MSTQW) found that the gravitational self force (in vacuum, and with no external forces) is due only to 
the tail part of the Green's function associated with the perturbation equation in the harmonic gauge. That is, the 
"instantaneous" part of this Green's function yields no contribution to the self force. [The formal expression for the 
gravitational self force as derived by MSTQW is given in Eq. (^ below.] 

We may point out here two aspects in which the gravitational self force differs from its electromagnetic (and scalar) 
counterparts. First, the gravitational self force contains no local term analogous to the electromagnetic (or scalar) ALD 
force (which is proportional to the time derivative of the particle's acceleration), as there is no acceleration associated 
with the purely gravitational field. A more fundamental difference concerns the very nature of the gravitational 
self- force as a gauge-dependent entity: When the metric is subjected to a gauge transformation (i.e., an infinitesimal 
coordinate transformation), additional terms emerge in the particle's equation of motion, which correspond to a change 
in the effective self force. (If fact, in the perturbative context the self force is a "pure gauge" entity, in the sense that 
it can be locally eliminated by a suitable choice of the gauge.) To provide any meaningful gauge-invariant physical 
information, the self force must therefore be supplemented by the metric perturbation to which it corresponds (which, 
of course, contains all information about the gauge). In this regard, an essential point is that MSTQW's prescription is 
formulated within the framework of the harmonic gauge. A priori, it is only in this gauge — the harmonic gauge — that 
the force is guaranteed to be well defined and finite. The analysis presented in the current paper remains within the 
framework of the harmonic gauge (namely, we shall discuss the "harmonic gauge self-force"). A systematic study of 
how the self force behaves under a change of gauge shall be presented in a forthcoming paper [Q. 

The formal expression obtained by MSTQW for the gravitational self force was sometimes considered impractical 
for actual calculations, as it was unclear how one should evaluate the nonlocal tail term in general cases. Also, to 
apply this expression, one encounters the problem of calculating the metric perturbation (and Green's function) in 
the harmonic gauge, for which perturbation formalism has not been fully developed as for other customarily used 
gauges (e.g., the Regge- Wheeler or radiative gauges). The first actual calculation of the gravitational self force based 
entirely on MSTQW was recently presented by Pfenning and Poisson who considered the motion of a mass 
particle (as well as scalar and electrically charged particles) moving in a weakly curved region of spacetime. In our 
present manuscript, the attempt is made to present a practical method for direct implementation of the MSTQW 
result in strong gravitational field. 

It should be commented that other approaches to the gravitational self force, not directly relying on the MSTQW 
analyses, were also taken recently by a number of authors. Lousto fll] proposed a scheme based on the zeta- 
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function regularization technique, to allow self force calculation in strong field (Lousto's scheme is similar to the 
method presented in this paper, in that both methods employ a multipole mode decomposition of the gravitational 
perturbation). Another method for extracting the finite part of the self force was proposed by Detweiler jllj. Most 
recently, Nakano and Sasaki (Tsj carried out a weak-field calculation of the self force in a Schwarzschild background 
by evaluating the tail part of a Green's function. It was assumed, however, that the correct force could be derived 
from a Green's function associated with a certain Klein-Gordon type wave operator, instead of the harmonic gauge- 
related Green's function as required in the MSTQW regularization scheme. It is unclear yet whether or not the 
above regularization schemes produce, in general, the results that would have been obtained by a strict application 
of MSTQW regularization. 

As we already mentioned, the calculation scheme to be presented in this manuscript is based on MSTQW's formal 
expression for the gravitational self force. It employs a technique — the mode-sum regularization scheme — previously 
introduced and tested for the scalar self force. In general, the application of this scheme consists of two 

essentially independent parts: in the first part, one expands the perturbation field of the particle into multipole modes 
(tensor harmonic modes in the gravitational case, e.g.), and derives the "bare" self-force Z-modes associated with the 
Z-modes of the metric perturbation. The Z-mode perturbation is finite and continuous even at the particle's location, 
and the corresponding Z-mode self-force maintains a finite value as well (although it usually suffers a discontinuity 
across the worldline). However, the sum over the bare force's /-modes turn out to be divergent. In the second part of 
the mode-sum scheme, certain regularization parameters are calculated analytically, by a local analysis of the Green's 
function at large I and small spacetime separations. These parameters are then used to regularize the divergent 
sum over bare force's Z-modes. This calculation scheme completely relies on MSTQW's results, and contains no 
further assumptions as to the appropriate regularization (although it does contain certain assumptions concerning the 
mathematical behaviour of various quantities involved in the analysis) . 

Previously, the above mode-sum scheme was fully implemented in several test cases: A scalar charge held static 
outside a Schwarzschild black hole A scalar charge in a circular orbit or one in radial motion 
Schwarzschild spacetime; and the motion of scalar or electric charges on the background of a massive shell (EOT 
cently, Burko and Liu first applied the mode-sum scheme for a static scalar charge in Kerr spacetime (however, an 
analytic foundation for the scheme has not yet been established in the Kerr case). So far, the mode-sum regularization 
scheme has never been applied for calculating the gravitational self force. 

The arrangement of this paper is as follows. We start, in Sec. ||, by summarizing and discussing MSTQW's 
result for the regularized gravitational self force, which is the basis for our calculation scheme. Specializing to the 




Schwarzschild geometry, in Sec. Ill we expand the (harmonic gauge) Green's function in tensor harmonics, and obtain 



a set of equations for its various multipole modes (i.e., a set which does not couple different I, m modes, although it 



does couple different tensor harmonic components). In Sec. IV we then introduce our mode-sum scheme as applied 
to the gravitational self-force problem. To demonstrate the applicability of this scheme, in Sec. ^ we implement it in 



two simple cases by explicitly calculating all necessary regularization parameters. Sec. VI contains a summary of our 



results and a discussion concerning the application of our scheme to other, more general cases. 

II. WORLDLINE INTEGRAL FORMULATION OF THE GRAVITATIONAL SELF FORCE 

We consider a pointlike particle of mass m moving in the exterior of a black hole with mass M 3> m.^ Let the 
metric 17^,^ describe the black hole background geometry, which we assume to be a solution of the vacuum Einstein 
equation (later we specialize our discussion to the Schwarzschild spacetimes). Let also z"(t) represent the particle's 
trajectory, and u°'{t) = dz^/dr stand for its four- velocity. The particle produces a small perturbation hap gafj to 
the background geometry. The Einstein field equation for the metric gap -\- hap, linearized in h^fj, takes the form 

Uha0{x) ~ 2R^afi''h^Ax) = -16^™ / {-g)-^^H^[x - z{T)]uo.{T)up{T)dT, (1) 



where we have introduced the "trace-reversed" metric perturbation 



ha0 — hap — -gaph (2) 



^The extent to which the concept of a pointhke particle makes sense in the context of the self- force problem is discussed in 
Ref. i. 
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(with h = g^'^h^y), and where we have set the harmonic gauge condition h^^ -^u = 0. In the perturbation equation 
(|l|), □ stands for the covariant D'Alambertian operator, Rap-yS is the Riemann tensor in the background geometry^, 
g stands for the determinant of the metric gap, and 5 is the Dirac delta function. The particle does not move along a 
geodesic of gapi as it interacts with its own field /ia/3- Phrased differently, to the perturbation hap there corresponds 
a self- force i^", in terms of which the particle's equation of motion is given by 

ma"(r) = F"(r), (3) 

where a" = u^^^u^ is the particle's four-acceleration, with a semicolon denoting covariant differentiation with respect 
to the background metric gap- 

Obviously, the perturbation hap diverges on the worldline itself, and the "bare" self force associated with this 
perturbation (as defined below) diverges as well. To obtain the physical equation of motion, one must appropriately 
regularize the self force. The analyses by MSTQW j6|Q| provide a formal expression for the regularized self force, to 
order 0(m?), in terms of a worldline integral over derivatives of a retarded Green's fimction. It is found that, in a 
vacuum background, F°' is solely due to the tail part of the Green's function: 

F''{T)=m'k''^'<' f Gp^p>^>,s[^^'{T)-z^'{T')]uP\T')u^\T')dT'. (4) 

Here, the "kinematic" tensor k"^^^ is given by 

I^cPjS ^ ^gO'S^Pyi _ g'^Pu^u^ _ + ^u°'g'^''u^ + ^g"''/^ (5) 

(understood to be evaluated at the particle's location), so designed to assure that the self- force has no component 
tangent to the worldline (i.e., F"ua = 0). This guarantees that the mass m maintains a constant value along the 
worldline. The quantity Gp-yp'-y' is the two-point retarded Green's function associated with the wave operator given 
in the left-hand side of Eq. (]^). It satisfies^ 

OGaPa'P'ix;x') - 2R^^ap''{x)Gi_,t,a'P'{x\x') = -167t{- gy^ ^"^ ga' (aix , x')g p)p, (x , x')5^ {x - x') = Zapa'p', (6) 

with the supplementary causality condition Gf^ua' p'{x]x') — whenever x lies outside the future light cone of x'. In 
this equation, the D'Alambertian is taken with respect to a;, parenthesized indices indicate symmetrization, and gaa' 
is the bi- vector of geodesic parallel displacement defined in Q . In what follows we shall need only the value of this 
bi-vector in the "coincidence" limit: \mix—,x' gaa' = gaa'- Note that the bi-tensor G ^i,a' p'{x]x') has, in general, 100 
independent components (compared with 16 components in the electromagnetic case, and only one in the scalar case). 
We also mention that the trace-reversed metric perturbation itself is constructed from the Green's function according 
to 

/OC 
GaPa'P'[x^]Z^^{T)]u"'{T)u''\T)dT. (7) 

To avoid confusion, it is worth commenting here about the different notation previously used by different authors 
in writing Eq. (^: The trace-reversed Green's function, denoted here by Gapa'p', is denoted G^®*^,^ by MST [||, 

^aba'b' B^Hl ' ^^'^ Qapa'p' by Pfenning and Poisson [|ll[. Note that Pfenning and Poisson use a different 

normalization for the Green's function (it contains an extra factor 4 with respect to either and the normalization 
used here). Also, note that in Ref. Q the self- force is expressed in terms of only the tail part of the Green's function, 
denoted therein by Vapa'p', with no contribution from its instantaneous part. This expression, however, is actually 
identical to the above Eq. (Eh, as the worldline integral involved is cutoff at r^. 



■^We use here the convention of Ref. ji^l for the Riemann tensor. Note the different convention used by Mino et al. in g. 

^ Note that although the wave operator defining the Green's function Gapa'p' indeed originates from the perturbation equation 
in the harmonic gauge, the Green's function itself does not satisfy the harmonic gauge condition, as one can directly verify. (To 
see this, note that the delta-function source does not satisfy the conservation law — a vanishing covariant divergence — as does 
the source for the metric perturbation itself. Consequently, the harmonic gauge condition is not consistent with the Green's 
function equation as it is with the perturbation equation.) Note that Eq. (2.12) in Q is therefore erroneous. 
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The Green's function G^,ya'p' plays a central role in our analysis. It is important to notice that, based on MSTQW 
analysis, Eq. (0) is guaranteed to yield the correct, "physical" force, only when using the Green's function defined 
through Eq. (^. The expression given in Eq. may fail to represent the physical self force, and may even yield an 
indefinite result, if a different Green's function is used. 

For future use, it is useful to write the Green's function equation (j6|) in the (non-covariant) form 



Here, Ds stands for the D'Alambertian operator acting on a scalar function. 



(8) 



(9) 



(where F^^ are the connection coefficients), and the ^„^'s are certain differential operators, of the first order at most, 
which describe how the various components of the Green's function couple to each other. In a given coordinate system 
we have 



A^"" G — 2n^ 



"2r('^(„GA/3),i/) 



G 



pA pp 



(10) 



where the source-point indices a'/?' have been suppressed for brevity. In the Appendix we give explicitly the operators 
At^a for the Schwarzschild background case (in Schwarzschild coordinates). 



III. TENSOR HARMONIC EXPANSION AND REDUCED EQUATIONS 

Accepting Eq. (jj) as the basic expression for the gravitational self force in curved spacetime, the main concern 
remains how to implement this expression in actual calculations. One may start by considering limiting cases, as the 
weak field or slow motion scenarios, in which Eq. (^) could be applied in a direct manner. This, indeed, was done 
recently by Pfenning and Poisson in Ref. ||ll|. However, in considering realistic black hole spacetimes, one ultimately 
wishes to apply Eq. (Q) for strong field self force calculations. Here, the introduction of a mode-sum scheme seems 
inevitable. A mode-sum decomposition approach is necessary for a dimensional reduction of the problem (in the 
usual manner), but it is especially beneficial in the context of the self force problem, because of the following reason: 
Whereas the metric perturbation diverges on the worldline, its individual modes, as well as the corresponding force 
modes, maintain a finite value even at the location of the particle. Thus, in exploring the behavior of the individual 
self force modes, one avoids dealing with divergent quantities. Still, introducing a mode-sum scheme for calculating 
the self force is not a straightforward task: the self force modes each carries a mixed imprint of both the tail and 
instantaneous parts of the worldline integral, which results in that the sum over modes usually turns out to diverge. 
A carefully designed scheme for regularization of the mode sum is thus necessary. The introduction of such a scheme 
is the main target of this paper. 

We first consider the multipole-mode decomposing of the Green's function. As in the rest of this paper, we focus 
on the spherically symmetric Schwarzschild black hole background, with a line element given by 

ds^ = -f{r)dt^ + f-\r)dr^ + r^{d0^ + sin^ Bdip^), (11) 

where /(r) = 1 — 2M/r and AI is the black hole's mass. Throughout this paper we use Schwarzschild coordinates 
t, r, 9, If, relativistic units (with G = c = 1), and metric signature — \ — \ — h. 

Any (sufficiently regular) symmetric covariant tensor of rank two. Tap, can be expanded on a 2-sphere in the form 

oo / 10 
;=0 m=-l 1=1 

where are scalar coefficients and Y^^^^"^ are the Regge- Wheeler- Zerilli tensor harmonics p^-p7[|. In Schwarzschild 

coordinates t, r, 9, tp, the set of tensor harmonics i^^^""™ is given by|^ 



^We adopt here the orthogonal set introduced by Zerilli [M, though we use a different notation for the basis tensors: The 
symbols Y^^J'"" . . . Y'j^°^'"' are used here instead of Zerilli's aj^, a^*^, a;™, b[^, him, gim, tim, c[^, cim, and dim, respectively. 
Note the sign error in Eq. (A2j) of Ref. H). 
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where y'"'(0, </?) are the scalar spherical harmonics, s = sinO, X = {I — + 2)/2, and the operators Di and D2 are 
given by 
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Di = 2(de - cot0)d^, D2 = dee - cot0 dg - s '^d^^. (15) 
The tcn0 tensor harmonies of Eqs. ( [l3| ) and (|lj) eonstitute an orthogonal set: 

dl]ry"^r;'^'^[r«""]*yij''™' = fc«%<5„,(5™„, (for j = 1, . . . , 10), (16) 



where 77"'^ = diag[— l,l,r ^, (rs) an asterisk denotes complex conjugation, and the integration is carried over a 
2-sphere (no summation over i is implied on the right-hand side). The coefficient fc^*) takes the value —1 for i = 2,4, 8, 

and +1 otherwise. The seven harmonics Y^^J^"^ . . . i^^^'™ constitute a basis for all symmetric covariant tensors of even 

parity, while the remaining three harmonics Y^^^"^ . . . y^'^^"*'™ span all odd parity tensors. Recall that the odd and 
even parity parts of a tensor (e.g., a metric perturbation) are uncoupled, and can be treated separately p4[ . Using 
the orthogonality relation ( p^ ) , one easily constructs the scalar coefficients of Eq. (^2[) through 

T«'™(r,t) fcW y rff7^"M^/3^[yW'™]*r„^. (17) 

Now, the Green's function Gaf3a'f3'ix,x') is a bi-tensor. It transforms like a tensor at the evaluation point x (when 
the coordinate transformation is carried out holding x' fixed), and it also transforms as a tensor at the source point 
x' (when the transformation is performed with fixed x). Regarding the Green's function, for a while, as a tensor at 
the evaluation point x, we may expand it in tensor spherical harmonics at that point, as in Eq. (p^). We write 

G^0^,^>ix,x') = irr'r' E gI^^^^C^, ^J^^U^'P'^ (18) 

where G'"^}j^ are the multipole expansion coefficients (independent of 6 and ip), the quantity G^^^,^, is the one 
resulting from formally summing over i and m, and the radial factor {rr')~^ is introduced for later convenience. 
We further need to expand the (bi-tensorial) source term of the Green's function equation (||) in tensor spherical 
harmonics. This expansion takes the form 

Z^p^^p- = irr')-' E ^S'"(^' ^i?""' (19) 

where, using Eq. ( p^ with Eq. (|6|), the coefficients are found to be given by 

- -167rfcW<5(t - mr - r0.gp.(„K^05^>'(^')/'^V')^^''''(^')[>;''i'™(^^')]* (20) 



(with fl' standing for 9',ip'). Note here that the coefficients zj^jj^T depend on the evaluation point x only through 
the delta functions]^ 

Eq. (||) is now separable into multipole modes l,mhy means of expansions ( [l^ ) and (|l9|): By substituting these 
expansions in Eq. (^) one obtains a set of equations which indeed couple between the ten functions G^^r^,'"^^^''™ (for 
given l,m), but not between the different multipoles l,m. To write the equations for the various multipole modes of 
the Green's function, it is convenient to introduce the Eddington-Finkelstein null coordinates 

v = t + r*, u = t-r* [where drVdr = /"^(r)], (21) 



^ Note that for / = and 1 = 1 there are actually fewer independent tensor harmonics: There are only 4 independent harmonics 
at Z = (the harmonics i = 1, 2, 3, and 6), and 8 independent harmonics for each of the three I = 1 modes (the harmonics i = 7 
and 10 vanish identically at / = 1 for any value of the azimuthal number m). 

® In fact, the naive construction of the coefficients Z^^i]^ yields an expression involving the hi- vector [as in Eq. The 

form (ji^) is then obtained by noticing that the coeflicients ^^V^™ transform like scalars at the evaluation point x. This allows 
us to prime all tensorial indices of the various factors involved in Eq. (^^ (i.e., take these factors to transform like tensors with 
respect to the source point a;'), without affecting the value of Zj^)'^ . The presence of the delta functions then further allows 
us to take all tensorial factors in Eq. (pol) to be functions of only the source point coordinates x' . 
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and the time-radial operator 

Di^d^., + V^{rl (22) 

where 

VUr)^Ul,'l±ll). ,23, 



4 \ r 

with / = /(r) and /' = df{r)/dr. The operator D\ is the familiar wave operator associated with the ^-mode of a 
massless scalar field in Schwarzschild spacetime [note the relation rUg[(f){r,t)Y^'^'^{^)/r\ = —Af^^{r)D\4){r,t)Y^"^{Vl), 



where (j)(r, t) is any function] . The equations for the various modes G|^*?i™ then take the form 



DiG^'P+T^G^j;' = Sf;- (sum over j), (24) 
where the source term is given by 

S^:^t' = -[/(0/4]^S™ - 8fc«^<5(z; - v')5{u - (25) 

and where are differential operators that couple between different i values. These coupling terms read (suppressing 
the indices a'/3' and Im for brevity) 



iSg(^) = i//'G« - \ {r + ^ff'/r) G(i) + I {ff^ - 2fr) G(3) _ iV2 //' [^Gf + r-'fG^'^) , (26a) 

lj-]G^'' = \ir - ff" + 2/V^')G(^) + \V2^fr (/-^G^ + Gf) l^lifTT) r'^ fG^^\ (26b) 

xgG(^) = -lffG[? + 1(8 - r + 4//'A)G^=^) + lif-'f" 2/-i/")G(i) - lV2^r'fGf 

-lV2r-^f^l{TTT)G^'^ + \V2{r/r - 2f/r')G^^\ (26c) 

igGf-'-) = Iff (g(,^) - 2r-iG(4) + ^Gf) - iyi(7TT) (/V^')G(^), (26d) 

i;;)G(^-) = -i//'G(:;) - 1/ (/'A - 4/A^) G(^) - i(/Vr^)v/^I(7TT) G(3) - \f-'f'^G[t^ 

+ fv/IOTT) G(6) - ^/2AG(^)l , (26e) 
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igGW = ^(//r^)(l - 2r/')G(«) - ^(/'A)G(i) - ^/^(2//r2 - /7r)G(3) + l,/l{rTT){f/r')&'\ (26f) 

^ (26g) 



xfJ^Gf-'") = i//' (G(«) - 2r-iG(8) - iGf^ (26h) 
X;2g(^^ = -Iff'Gl? - If {f'/r if/r') G^^) + \f-'f'^Gf - ^if y)V2XG(''\ (26i) 

^(10)^0) ^ _1(^/,2)^(10) _ l(;2/^2)^^(9)_ (26j) 

The separated equations ( p^ ) have the convenient property that no coupling between the various modes occur in 
the main parts of the equations (i.e., the parts containing second derivatives). Coupling between different i-modes 
comes into action only through the terms, which contain one or no derivatives. Note that the even parity modes 
(z = 1 ... 7) do not couple at all to the odd parity modes {i = 8 . . . 10). The two mode types form a disjoint set of 
equations, as one would expect. 
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IV. MODE SUM REGULARIZATION 



Without loss of generality, let us take the point along the particle's trajectory where we wish to calculate the self 
force to be at T = 0. Based on Eq. (Q), we may express the self force as 

F„ =i^^barc) _^anBt)^ (27) 

where 

.0+ 

F^bare) ^ ^2j^J^5 / Q ,s[z>' {t = Q)' Z^' {t')]uI'' {t')v?' {T')dT' = mkj^%^,s (28) 

J —OQ 

is the "bare" force associated with the metric perturbation hap [the second equality here stems from Eq. (^], and 



^(inst) ^ Ymi SFj;'^ = hm 

£-.0+ e^0+ 



+ 



m2fc„/37<5 / Gf3-,f3>Y.5[z^'iT = 0);z^'{T')]ul^' {t')u-'' {T')dT' 



(29) 



is the "instantaneous" part of the force. The quantities F^'^'^^'^ and i^^™''*^ — both involving integration through the 
particle's location — are singular and so poorly defined as they stand. For definiteness, we may redefine the integrands 
appearing in Eqs. (|2^) and ( |29| ) as vector fields in the neighborhood of the evaluation point, and later be interested 

in their behavior on the worldline. Note, however, that the difference Fj^'^'^'^'' — Fq'"''*-' does yield a definite finite value 
at the particle's location: According to the analysis by MSTQW, this value represents the physical self- force Fa. 

To introduce our mode-sum regularization scheme, let us denote by pj^^'^'^'^") and SFa ^\ respectively, the contri- 
butions to F^^'^'^'^^ and SFa ^ coming from the ^-mode of the Green's function [these two quantities are obtained by 
replacing the Green's function Gpjp'^' in Eqs. (|28| ) and ( |29| ) with its /-mode G^^^,^, defined in Eq. We may 

then express the self force as a sum over /-modes, in the formQ 

Fa = lim^ J2 (i^o^*"""^ - SFi''>'^ . (30) 

An essential feature of our scheme arises from the fact that the /-mode of the metric perturbation haf3 is everywhere 
finite and continuous; it remains finite and continuous even at the location of the particle, where the overall per- 
turbation diverges. Consequently, the /-modes of the bare and instantaneous forces, pj^^'^'^'^") and SFa^\ turn out to 
attain finite values. This behavior has been analyzed and demonstrated in the analogous scalar self force problem 
| p5| , p^ , and is equally valid in the gravitational case as well (we demonstrate this in the next section) . As in the scalar 
self-force model, the two quantities p^^'^'^'^') and 6Fa^^ are discontinuous through the particle's location [regarding the 
integrands appearing in Eqs. ( |2^ ) and ( p9| ) as vector fields in the neighborhood of the particle]. That is, each of these 
two quantities attains different (finite) values if calculated by taking the limit r {z'^)~^ or else the limit r {z'^)~ ■ 
[Later we assign to pj^^^'^^') and SFa'^'' the labels -I- or — to indicate weather they were calculated from r [z'')'^ 
or rather from r — > {z^)~ , respectively.] Note, however, that the difference pj^^^'^"'' _ SFa ^'' [as well as the sum over 
modes in Eq. (pO|), producing the regularized force Fa] does not depend on the direction from which the limit is taken. 

Although each of the bare modes pj^^'^'^'^') yields a finite contribution to the self force, the infinite sum over pj^^^'^'^^ 
diverges, in general. This is easily demonstrated already in the simple case of a static mass in fiat space, a case studied 
in the following section (see also [^ for a discussion of this point in the analogous scalar case). To carry out the 

regularization procedure, one seeks an (e-independent) function H^J^, such that the series J^ii-^a ^^''''^^ ~ H'^) would 
converge. Once such a function is found, Eq. (^0|) can be written as 

oo 

= E - ffi^) - Dt, (31) 

1=0 



It is assumed here that the differentiation and the integration involved in constructing F^^'^'^'^^ and SFa ^ out of G/3^p>jr can 
both be performed term by term with respect to the sum over /. 
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where 



D^^ lim V - . (32) 



. .0+ • 

The parameter is well defined, since, by construction of H^^ , the difference 5Fa'^^^ — yields a convergent 
sum over I [note that the two quantities p^^^'^^°^^ ^nd 5Fa^^^ must bare the same singular behavior at large Z, as 
their difference yields a convergent sum over I — see Eq. (|30|)]. Also note here that the limit e — > and the sum over I 

should not be interchanged; otherwise, the crucial contribution from 6Fa^''^ would be lost. 

Eq. constitutes the basic expression for the gravitational self force in the framework of the mode sum regular- 
ization approach. The implementation of this expression for calculating the self force at a certain point along a given 
trajectory involves two (essentially independent) parts: In the first part, one should first calculate (numerically, in 
general) the l-modes of the bare metric perturbation (in the harmonic gauge) and then use these modes to construct 

the bare force modes Fa^'^'^''\ In the second part of the calculation procedure one should obtain the "regularization 
function" i?^. In principle, this function should be constructed by exploring the asymptotic behavior of the bare 
modes as / ^ oo. It is more convenient, however, to read this large I asymptotic behavior from the quantity 
5Fa'^\ which is strictly local (recall that p^^^'"^^ ^nd 6Fa^^ bare the same singular behavior at large I). To this end, 

one merely needs the asymptotic behavior of SFa''^ in the immediate neighborhood of e = 0. This allows one to derive 
i?^ (and later also Da) using local analytic methods, as we shall demonstrate in the next section. 

In general, the (one-sided values of the) bare modes p^^'^'"^^ and SFa'^'' are found to diverge at large Z as oc L This 
was demonstrated for various trajectories in the scalar self force problem, and will be demonstrated below in the 
gravitational case as well. To construct the regularization function iJ^ so as to regularize the mode-sum in Eq. (pi]), 
one should therefore derive the three leading-order terms in the l/l expansion of 5Fa ^\ It appears more convenient 
to use an expansion in powers of the new variable 

i = / + l/2. (33) 

Denoting the coefficients of the leading-order terms by Aa, i?Q, and Ca, we shall have, in general, 

Hl=AaL + Ba+Ca/L. (34) 

The implementation of the mode sum regularization scheme therefore amounts to (i) Calculating the bare modes 
^^(baro) ^|^]_^^g usually douc numerically); (ii) Deriving the four "regularization parameters" Aa, Ba, Ca, and Da 



(by local analytic methods); and (iii) summing over I using Eq. ( |3l| ) to obtain the regularized force Fa- 

The scheme described here is based on the prescription (^, which is formulated within the harmonic gauge. It 
therefore requires that the bare force modes be obtained from the metric perturbation in the harmonic gauge. This 
poses a problem from a practical point of view, as a separable formalism for the metric perturbation in the harmonic 
gauge has not been well developed as it has for other gauges [especially the Regge- Wheeler (RW) gauge] We shall 
deal with this gauge problem in a forthcoming paper [nO| , where an attempt is made to rephrase the scheme in terms 
of other gauges, such as the RW gauge. 

Before proceeding to discuss some simple applications of the method proposed here, we should comment on a 
certain technical subtlety which we avoided so far in our discussion. As we thoroughly discuss in Rcf. jl^ (in the 
context of the scalar self force problem), the Green's function does not admit a convergent multipole expansion, as 
a result of its being singular along the light cone of the source point. As a consequence, the modes jt'^^'^''™' (and 
SFa^^) contain certain terms which oscillate rapidly at large /, rendering the sums over I in Eqs. (pO[), (pl|), and 
(^H) non-convergent. In Ref. we justified throwing away these divergent oscillatory terms. To formalize the 
omission of these terms, we therein introduced a new summation and limit operations (the "tilde-summation" and 
"tilde- limit" ) which eliminate any oscillatory divergent terms while preserving the monotonic behaviour. The same 
problem — with the same solution — persists in our current gravitational case. However, to avoid complexity in our 
current presentation, we shall not attempt to indicate explicitly where a tilde-operation is to be applied [as in above 
Eqs. (|30|), (^, and (^]. In the analysis to follow, a tilde-summation or a tilde-hmit will be imphcitly used when 
appropriate. 



This problem becomes more acute when dealing with Kerr spacetime, for which a separation formalism for the metric 
perturbation has been developed so far only in the radiation gauge |28|. 
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V. SIMPLE APPLICATIONS 



We now demonstrate the apphcabihty of the above calculation scheme in two simple test cases. First, we consider 
the trivial case of a static point mass in flat spacetime. This would provide a simple test case (where the result is 
obvious: a vanishing self force) against which we may check the validity of our scheme. We then move on to the 
Schwarzschild spacetime, and consider a freely falling particle on a radial geodesic. In both cases we construct all four 
necessary regularization parameters. For simplicity, when considering the second case we shall focus on calculating 
the force at a presumed turning point of the geodesic (i.e., where dz^ / dr — 0), for which case the calculation becomes 
considerably more simple (see below). We emphasize that our calculation and results apply equally well for either 
weak or strong field. The a pph cation of the scheme in more realistic cases (ones of greater astrophysical relevance) 
will be presented elsewhere pyj. 



A. Static particle in flat space 



We consider a static particle of mass m in Minkowski spacetime. The particle is located at r = tq, in a certain 
spherical coordinate system r, 0, (/?. (We adopt here spherical coordinates in order to make the calculation more 
closely related to the Schwarzschild case discussed later.) The four- velocity of this static particle is m"(t) = 6f at 
all T, and thus only the Green's function components Gp^t't' take part in constructing the force through Eq. (jj). 
Considering the system of equations (^) for the 10 functions Gpjt't', one finds that the source Za/Bt't' is nonvanishing 
only for a/? = tt. [To see this, note in Eq. ( pO| ) that for a'/3' — t't' the only contribution comes at i = 1. Then, the 

only contribution to 5^^^'"™ in Eq. ( |l9| ) is at af3 = tt.] Since all interaction terms A^^f^Gut't' appearing in Eqs. (||) 
vanish in Minkowski spacetime [see Eqs. (Al)-(AlC), with /' = 0], we find that only the sourced component Gttt't' 



takes a nonzero value, while all other components Gp-yt't' (which satisfy homogeneous equations) vanish. Thus, for 
a'p' = t't' , the system of equations (H) reduces to a single equation for the quantity Gttt't' = G: 

□,G = Zttt't' = -Wni-gy^^^S^x - x'). (35) 

Note that in the simple case considered here — a static particle in flat spacetime — only one Green's function component 
out of 100 actually takes part in the computation of the self force. 

The self force can now be constructed from Eq. (||) by setting u°' = (5f and u" = <5" , and recalling that all 
component Gp^t't' but Gttt't' vanish. We find (evaluating the force at r = without loss of generality) 

1 /■°" - 

F^^-m^ G,„[z^(t = 0);z^(T')]dT' iova = r,9,ip, (36) 

^ J — oo 

as well as Ft = 0. [For a — t, the integrand in Eq. (^) in identically zero. This is a trivial result for a static particle, 
as the force is known to satisfy the normalization condition, Fau" = 0.] 

The Green's function equation (|3^) is exactly the same as the Green's function equation for a scalar field [compare 
with Eq. (4) in Ref. [p^], apart from a relative numerical factor of 4 on the right-hand side (the source in the 
gravitational case is greater by 4 than the source in the scalar case). In addition, the construction of the self force 
from the Green's function through Eq. ( p6| ) is exactly the same as the construction of the tail part of the force in the 
scalar case [compare with Eqs. (12) and (13) of Ref. |l^], apart from a relative factor of 1/4 on the right-hand side, 
which compensates for the extra factor of 4 in the Green's function equation. We may conclude, in particular, that 
the bare modes of the self force acting on our static particle of mass m are equal to the bare modes of the scalar self 
force acting on a static particle of scalar charge q = m. It is then possible to simply use the results already obtained 
in the scalar case: For the regularization parameters we obtain (taking the Minkowski limit of the results described 
in Refs. ll^,|lH| and replacing q ^ m), 
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^J = T^5^, Ba.^--^S'^, Cc^^Do^^O. (37) 
As the average of the two sided values of Aa vanishes, we finally obtain from Eq. ( pl| ) 

oo 

= E (i^a^''""^ - i?a) , (38) 
1=0 

where ^"Hbarc) ^ J^!(baro)+ _|_ ^;(baro)-j^2. 
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Now, in our current trivial case, the averaged bare- force modes appearing in Eq. (38) are easily calculated [ p^ : 
Solving first for the metric perturbation modesy, one finds that all components of the trace-reversed modes /i^^ vanish, 
except /ijj, which is given by — mr^'-^^rQ (for r > ro) and = mr'-r^^^^ (for r < tq). Then, using the second 
equality of Eq. (p8|), one obtains ~ —{I + l)77i^rQ ^ and — lm?rQ^ , yielding — —m'^/{2rl) (with all other 
components vanishing). Thus, the averaged bare modes of the force are found to be /-independent, each identically 
equal to the regularization parameter Ba ■ Consequently, one finds that each of the terms in the sum over I in Eq. 
( p8| ) vanishes independently, with an obvious (and expected) vanishing of the overall self force. 

B. Radial geodesic motion in Schwarzschild, at a turning point 

Let us now consider a particle moving in a radial geodesic in Schwarzschild spacetime. For this orbit we have, 
identically, = u'^ = 0, and therefore only the Green's function components Ga/st't', Gapt'r' = Ga/Sr't', and Ga/Sr'r' 
would take part in constructing the force through Eq. (Q). Now, for any given combination of a'/3', Eq. (H) constitutes 
a set of 10 coupled equations for the 10 independent quantities Gapa'p'- Considering the three sets of equations (|^) 
with a' [3' = t't\ t'r' , and r'r', we find that the source Zapa'p' is nonvanishing only for a[3 = a'(3' (e.g., in the set 
of equations for Gafit't', only the equation for Gttt't' is sourced). Similarly, Eq. ( |24| ) forms, for any specific value of 
a'/?', a set of ten equations for the ten tensor-harmonic modes G'^'7J3^■■^*'''''". The source S^^^^ for these equations is 
nonvanishing only at z = 1 for a' (3' = t't' , at z = 2 for a' (3' = t'r',r't', or at i — 3 for a'P' = r'r'. Clearly, the three 
odd-parity modes i ~ 8,9, 10, which are not sourced and also do not couple to any of the even parity modes z = 1, 2, 3, 
would all vanish (this is expected, of course, by virtue of our physical setup, which only includes an even-parity 
source). However, the four even-parity modes z = 4, 5, 6, 7, although not sourced in Eqs. (p^), do couple to the modes 
z = 1, 2, 3 and will therefore not vanish, in general. 

In conclusion of the above discussion, we find that in the scenario considered here — a particle moving radially on 
a spherically symmetric background — one has to deal with 3 sets of 7 coupled equations for a total of 21 nontrivial 
independent components gI^?^?': One of these sets (corresponding to a'f3' = t't') contains a source term only at z = 1, 
the other set (corresponding to a'P' = t'r' or r't') is sourced only at i — 2, and the third (for a' (3' = r'r') only at 
z = 3. To write these three sets of equations in a simple form, it is convenient to re-define our spherical coordinate 
system, such that the radial trajectory would be directed along the polar axis. In this spherical system, the Green's 
function (now sourced only at 9' — 0) would contain only axially-symmetric, m — modes. We then also introduce 
the new variables G|j?^,, defined (for a' f3' — t't' , t'r' , r't' , or r'r') through 

G% = 87Tgp,^^,{x')g(3,)^,{x')v'''^' V^'-"' a^'^Y'{e')G% (no summation over i), (39) 

where Y^{9) = y'^™=o and a'*-' = 1 for all i, except a'^^^ — i/V2. Then, each of the above three sets of equations (for 
a' (3' = t't' , t'r' , or r'r') takes the form 

DiGl^l,, +X«/g|^?; = q%6(u - u')Siv - v'), (40) 

where 

sts*;,, z-i, 

'^^^^ <! 5"?')' (41) 
0, z = 4,...,7, 

and j'^j^^ — for all except l|^^2)' = '^^'^^-^{2) ^'^'^ ^(j^2)' ~ '^(j)^'/'^''^''' Fi^^^Uy? to express the /-mode 

Green's function G^^^,^, in terms of the new variables G^*?^,, we substitute for G^^^^, in Eq. (p^), and consider only 
the Green's function components with aj3 — tt, tr, rr (we shall need only these three components in our following 
analysis). Recalling Y\e) = [{21 + 1) /{in)]^^^ Pi {cos 9) = [L / {2n)]^/^ Pi {cos 0) (where Pi is the Legendre polynomial) 



® Interestingly, in our case — a static particle in flat spacetime — the metric perturbation in the harmonic gauge exactly coincides 
with the one in the RW gauge. 
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and Pi{cos9') = 1 for 6' = 0, we then obtain for the three relevant components of the /-mode Green's function, 
evaluated at the polar axis (9 = 0), 

G\,^,p, = ^{rrr'f{r')LG%,, Gl^,p, = -2{rrr'LG%, = i{rr')-' r\r')LG%,. (42) 

[We may already, at this stage, evaluate the Green's function at the polar axis (where the particle is located), as the 
following construction of the r component of the force does not involve differentiation of the Green's function with 
respect to neither 9 nor ip — see Eq. (|3|) below.] 

To implement our mode sum prescription, we would now like to derive the regularization parameters A^, B^, Ca and 
Da ■ This task turns out to be considerably simpler (though not trivial) when the self force is evaluated at a presumed 
turning point of the radial geodesic. As a first non-trivial demonstration of applying our mode-sum regularization 
scheme, we hereafter focus on this special case. That is, we assume that the geodesic particle is momentarily at 
rest (i.e., m"" = 0) at r = ro, and calculate the local self force at that point. For simplicity, we shall consider in our 
calculation only the r component of the forcej^ Extension of this calculation to an arbitrary point of a radial geodesic 
is straightforward though rather laborious, and will be treated elsewhere |2^. On the technical side, the calculation 
of the regularization parameters very much resembles the calculation in the scalar case, which we describe in much 
detail in Ref. |]T6| . For this reason, we avoid here many of the technicalities involved in deriving the parameters, and 
refer the reader to Ref. for more details. In what follows we only give a very general description of how the 
analysis proceeds in our current gravitational case. 



As explained in Sec. IV, the mode-sum scheme's regularization parameters can be derived by exploring the behavior 



of the quantity SFa^^ [see Eq. (|29|)] at small e, using a perturbation analysis of the Green's function Z-modes at large-Z 
and small spacetime deviations. In our current problem of a radially moving particle at a turning point, the radial 
component SFr'^'^'' is constructed from the various components of the Z-mode Green's function by 



(with / and /' = df /dr understood to be evaluated at tq), which is obtained from Eq. (29) by setting — u = uf ^ 
and u* = f~^^^- 

To analyze the large I behavior of the /-modes G^^^,^, appearing in Eq. (^), we first introduce the "neutralized" 
spacetime deviation variables^ 

Ar = foL{r^ - r*o), 
Ar' = foL{r[ - r*o), 

At = foL{t-to), 
At' = foHto - t'), 

z^foL[iu^u')iv-v')f\ (44) 

where (r, t) is an off-worldline point in the neighborhood of ro, (r', t') is a worldline point in the past neighborhood of 
To, and /o = /^^^(7'o)/^o- We then consider the /-mode Green's function as being dependent only on L and the above 
"neutral" variables, and formally expand the quantities gI^}^, in powers of while holding the "neutral" variables 
fixed: 

oo 



To derive all necessary regularization parameters (including Da), it is sufficient to obtain an expression for the three 
leading-order terms in the 1/L expansion of the integrand in Eq. ([43|). Higher order terms in this expansion do not 



^ The angular components Fe and Fip are expected to vanish due to the symmetry of our physical setup, although we shall not 
attempt here to verify that our scheme indeed leads to this obvious result. At a turning point we also have Ft = 0, stemming, 
in a trivial manner, from the orthogonality relation Fau" = 0. 

^^The variables A^, A^/, At, and z are "neutral" in the sense that they each consist of a small, 0(e) spacetime deviation, 
multiplied by the large quantity L. The motivation for introducing this kind of variables stems from the detailed discussion in 



13 



affect the values of the regularization parameters . By analyzing the Green's function equations (^0|) , we now show 
that the contribution to SF^^^ conies at relevant order only from the three components G^jj/^/, Gj^j,j, , and G\^^,^,. In 
particular, it is shown that the terms of Eq. (^3|) involving Gj,^^/^, contribute (for any a (3') only at irrelevant high 
order. 

We start by substituting expansion (p5| ) into the set of Green's function equations (pol), and pointing out a few 
useful "rules of thumb"^: (i) When the operator acts on a mode G^*?^, it "lowers" its order by 2; namely, if 
G^'?^, oc 0{L") (for some n), then DgG^*,^^, cx 0{L"^^). (ii) A ^-derivative always acts to lower the order of G^?^/ by 
one: if gI^}^, oc 0{L"), then G^^jj^, ^ oc 0(L"+^). (iii) An r-derivative may lower the order by one at most when acting 
on a function of both z and A^, but does not affect the order when acting on a function of z alone.|^ (iv) The source 



terms oc S{u — u')S{v — v') appearing in Eqs. (40) are of order L^. (v) An immediate consequence of all above rules is 



now apparent from examining Eqs. (40): For any given z, a A: = contribution to a mode G^Vi, can only come from 



the source term (when non- vanishing) , whereas the coupling terms I^^^'G^?^, contribute only at higher (fc > 0) order. 

Consider first Eqs. ( ^0| ) for a' [3' — t't'. In this case, the only sourced equation is the one for g]}^! ■ After substituting 
expansion (^), we find by the above rule (v) that at the k = order this equation takes the simple form DlG^^j^j^^^^ — 
d{u — u')d{v — v'), with no contribution from the coupling terms Tjjj''G['ij! at this order. The solution to this equation 



is of order L° [see Eq. (49a) below], thus g[}J! oc 0{L°). Now, the mode G^?^/ is coupled to g[}^! through the g[1^ 



term in Eq. ( ^6b| ). By the above rules (i) and (ii) we thus find g|?j/ cc 0{L ^). On the other hand, the mode G[fJ! is 



excited only aX k — 2, through coupling with G^pji and g|?j/( — see Eq. (26c). We find further that G).?^^''' oc 0{L ^), 



=.(4,5); ^ nfT-3\ 



G|f(/ oc 0{L ^), and &J}! oc 0{L "'). It can be easily checked now that, up to 0{L ^), the first two equations of the 



set (24) — the ones with i = 1 and i — 2 — form a closed set of equations, with coupling to other modes affecting only 



at higher orders: 

Using Eq. ( p2| ) we may now evaluate the order of the /-modes G^^(,j, , and proceed to evaluate the order of the five 

integrand terms involved in constructing SF^^^ through Eq. (^). Using rules (ii) and (iii) and recalling u* oc 0(L'^), 
we find that the first integrand term in Eq. ( ^3[ ) is of order L^, the second and third terms are of order L^, and the last 
two terms (the ones involving G^.^^/^,) are only of order L~^. Since only integrand terms up to 0{L^) are necessary 
for calculating the regularization parameters, we conclude that the component G j^fif , would be of no relevance for 
this calculation. Up to the necessary order, the set of Green's function equations ( pjj ) therefore reduces, in the case 
a'P' = t't', to a closed-form set for g|,^'/ and g[^^J , given by Eq. (g|). 

Consider next the set of equations ( j40| ) for a'P' — t'r' . Here, only the i — 2 component is sourced, and we find 
G|?^f oc 0(LO). The modes i = 1 and i = 3 are sourced by G|?^f j, leading to g[}^J^'' oc 0(L ^). One similarly finds 
that the mode i = 4 is of order L^^, the mode i = 5 is of order L~^, and the modes i = 6,7 are only of order L^'^. 



Again, we use Eq. (|4^) to evaluate the order of each of the five integrand terms appearing in Eq. ( |43[ ) , this time for 
a'f3' — t'r' . We now recall, however, that u'' vanishes at tq, and is therefore of order m'' oc r/L. It is then easily 
shown that the only relevant integrand term is G^^^/j/ ^u* w'' oc 0{L^), while each of the four other terms contribute to 
the integral in Eq. (Bs) only at order L~'^, which is irrelevant for calculating the regularization parameters. Moreover, 

~ (2)1 

one finds that at the relevant order, the mode GJ,^, satisfies a single closed-form equation: 



In the following discussion we use a terminology according to which the "order" of a mode G^'?^, is determined by its 
expansion through Eq. (^), where the "neutral" variables are held fixed. That is, the "order" of G^V^, is L~^° , where fco is 
the smallest index k for which G)^;^,jj,] is non- vanishing. 

This behavior (iii) is special to a turning point. To see this, note that at r = ro we have dz/dr — (fo/ f)LAr' /z. Expanding 
A,,/ and z in powers of r about r(ro) = and defining the "neutral" proper time f as in Eq. ( |5C| ) below, we find, at f = 0, 
Ar' — f/tT^'^^^/^ ^11^ z = r (to leading order in 1/L). Consequently, if f{z) is some function of a certain order in 1/L, then 
the r-derivative df{z)/dr = [df{z)/dz]{dz/dr) = y~^[df{z)/dz] rr remains of the same order. 
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DsG^^ll = 5{v - v')5{u - u'). (47) 

(This equation happens to coincide with the scalar Green's function equation.) 

Finally, let us consider Eqs. ( ^ ) for a' j3' — r'r' . The only sourced mode is now cfr'^,, which is therefore of order 
L". However, recalling that at a turning point we have {u^ )^ oc 0(L~^), we find that the contribution from this mode 
to the integral in Eq. (^3|) is only at order L~^, and that the contribution from G')!?^, and G^,?^! is at still higher order. 

Therefore, no relevant contributions to SF^^^ arise for a'/3' = r'r' . 

In conclusion, when restricting our analysis to the case of a turning point, the problem of calculating the self 
force via our mode-sum scheme simplifies considerably: Instead of the three sets of seven equations each for the 21 
components required for an arbitrary point of a radial geodesic, one now has to deal with only the three equations 
( ^ ) an d (|47| ) (of which two are coupled and one is closed) for the three components g|/^(/, G['f^!, and C?!?^/'. 



Eqs. ( [461) and (|47|) are solvable in a perturbative manner, using the technique described in detail in Sec. V of Rcf. 
[^|. To apply this technique, one first writes 

G%=G'^X,eiu-u')e{v-v') (^ = 1,2,3), (48) 

where & is the standard step function, acting to confine the support of the Green's function to within the future 
light cone of the source point x' . Substituting this expression into Eqs. ( |4^ ) and (^^, one then finds that the new 

quantities G[^}j^, (treated as functions of v,u with a fixed source point v',u') must satisfy the homogeneous part of 

these equations anywhere at w > t;' and u > u' . One also finds that the value of G')^,^^, along the null rays v — v' 



and u = u' is uniquely determined from Eqs. (46) and (^7|). [Obtaining these "initial values" for the quantities G^V^/ 
involves the solution of a set of ordinary differential equations along v = v' and u = u' p9t . In the scalar case, the 
Green's function was found, in this way, to admit a constant value (of unity) along these initial rays [^ . In the 
gravitational case, the "initial data" are a bit more complicated, and will be given elsewhere |2^].] Thus, in effect, 
the above procedure converts the original set of Green's function equations into a characteristic initial-data problem 
for the quantities GJ^,, with a unique solution. This unique solutions reads, to the relevant order, 

+ { [{Ar - A.OK/I + /4)A. + /4A,,] + /|A2/4] Jo(z) + [f!{Ar + A,,0' - h] z^J2{z)/2 + f!z^Mz)/6 
+ [-/i/2(A2 - A^,) - /3(A2, + A,A,, + A^,) + /g] z Ji(z)} L-2 + 0(1-^), (49a) 

4'*/ = -r'f2{At + At,)Jo{z)L-^ + [hzJi{z) + frMz)] Ar'AtL'^ + 0(1-^), (49b) 

G[?^f = Jo{z) - /iz Ji(z)(A, + Ar,)L-^ + 0{L-^), (49c) 

where J„(z) are the Bessel functions of the first kind, of order n, and the /„'s are radial factors given by /i = 
i/-i/2(^/'_2/), h - r/'r^/^ h = r2[/-i(/')2 + /"]/12 + (/-r/')/2, h - r^f"/2, h = U + fi /2 + rf /2 + 1/8, 
h = 7-^2(2/1/2 _ j2)/8, and h = r^f-^f" - 2f-^f'f]/4 (aU evaluated at r = ro). 

The analysis now proceeds as follows: (i) We substitute the solutions (^9|) in Eq. ( ^2| ) to obtain the relevant /-modes 
^L/3q'/3'i ^^'^ then substitute these /-modes into Eq. (|4^). (ii) To be able to carry out the integration in Eq. (^3|), we 
next expand all a;'-dependent quantities now appearing in the integrand in powers of r about the evaluation point 
T(ro) ~ 0. (This procedure is described in detail, as apphed to the analogous scalar case, in Sec. VI of Ref. fl^.) (iii) 
We introduce the "neutral" proper time variable^ 

f = -{L/ro)T, (50) 

and write the integrand as an expansion in powers of l/L, with f held fixed. 
Following these manipulations, Eq. (E3|) takes the form 



Note the different notation used in Ref. ||l6|, where the "neutral" proper time variable r has been denoted by A. 
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rLe/ro ^ , 

dFy^'^=m' LHi"^^if) + Hi'\f) + Hj,^\f)/L + 0{L-^) df, (51) 

Jo L J 

where Hr^^ are certain functions of only f and ro (but do not depend on / otherwise). These functions all have the 
form of a sum over a few terms oc f*'V„(f), where fc,n € N [see, in the analogous scalar case, Eqs. (88)-(93) of Ref. 
1^ ]. The function Hr°^ has two different values, denoted in Eq. (|5^) by H^^^^ and , according to whether the 

derivatives involved in constructing Hr^'^ are taken from r or r , respectively. This kind of discontinuity, 

which shows up only at the leading order in the l/L expansion, results from differentiating the Q{u — u')Q{v — v') 
factor appearing in Eq. (|4^): the contribution coming from the light cone [through the S{u — u') or S{v — v') factors] 
depends, in its overall sign, upon the direction through which this derivative is taken. This effect is further discussed 
and illustrated in Sec. IVc of Ref. |l^ (see especially Fig. 1 therein). 

In terms of the functions Hi^^ , the regularization parameters are constructed by [ [T^ 



/•oo />oo /-oc 

Af^m^ iJ^o)±(f)df, Br=m^ i/W(f)df, Cr = HP{f)df, 
Jo Jo Jo 



(52) 



and 



Dr = -m^ lim y / LHI°^^ + H^^^ + iJ^^) (53) 

[Both functions iJ^"^^ and h'^'^^ can be shown |l^ to yield the same contribution to the integral in Eq. (|53|), which 
is why no ± sign has been assigned to the parameter Dr .] The evaluation of the integrals over f [and of the sum over 
I in Eq. (|3|)] is done in a manner completely analogous to the scalar case, as described in Sec. VII of |l^]. Here we 
merely give the results of this calculation: The parameters , Br, Cr, and Dr are found in our case — a mass particle 
at a turning point of a radial geodesic — to be given by 

At = T^(l-—] Br = -^, Cr = Dr = 0. (54) 



ro / 2r| 







We comment that the vanishing of the parameter Cr is necessary to assure consistency of our entire scheme: otherwise, 
the parameter Dr would have been indefinite (this point is explained in Sec. VIID of [p^).p| 

The values derived here for the regularization parameters find support from a recent numerical analysis carried 
out by Lousto |12 2^: Lousto has calculated numerically the bare modes of the force, Fq'^'^™', for a radially moving 



particle in Schwarzschild spacetime (as part of the implementation of a different regularization scheme), and found 
that these modes indeed show a large-Z behavior of the form indicated in Eq. (|3^) above. Furthermore, the analytic 
expressions derived here for the coefficients A^, Br, and Cr show a perfect agreement with the numerically-derived 
coefficients p9| . In addition, our result of a vanishing parameter Da confirms, in the cases studied here, Lousto's 
result based on his proposed zeta-function regularization. (Although Lousto's numerical calculations were carried out 
in a different gauge — in the RW gauge rather than in the harmonic gauge — this should not alter the values of the 
regularization parameters, as we explain in Ref. pO|.) Finally, we mention that the expressions derived here for A^ 
and Br also agree with the analytic values obtained recently using a different formalism, in which the self force is 
derived from the so-called Moncrief waveform [|2^ 

In conclusion of this section, we have found that the r component of the regularized self force at a turning point 
of a radial geodesic in Schwarzschild spacetime is given by above Eq. (|38|), where the parameter Br is given in Eq. 
(p^), and where p^^^'^'-'^'> qj-q the (sided average) /-modes of the bare force. These bare modes are derived [through 
the second equality of Eq. (p8|)] from the metric perturbation in the harmonic gauge. For practical use, it would be 
desirable to express our result in a more useful gauge (e.g., in the RW gauge). This shall be done in a forthcoming 
paper |l^, as part of a more general discussion of the gauge issue in the context of the gravitational self force problem. 



Interestingly, the parameter values (^4|) coincide with the values obtained for the scalar self force acting on a particle of 
scalar charge g = m, at a turning point of a radial geodesic in Schwarzschild spacetime — see Ref. jl6|. 

So far, agreement has been achieved for both one-sided values of Ar and Br, as well as for the average value of Cr- Currently, 
we obtain, using Moncrief 's formalism, that although the averaged Cr vanishes (as in our present analysis), the one-sided values 
of Cr fail to vanish. It is most likely that this preliminary result is erroneous. This point awaits further examination. 
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VI. SUMMARY AND DISCUSSION 



In this manuscript we have generahzed the scheme of mode-sum regularization, previously apphed only in the 
scalar case, to the problem of calculating the gravitational self-force on a mass particle. The proposed scheme offers 
a practical way for implementing the formal prescription developed by MSTQW, even in strong-field calculations. 

Within the mode-sum scheme, the basic formula for constructing the (harmonic-gauge related) gravitational self- 
force is given by 



= E ([J^a^'^"'"^]" - Aa,L - B„ - C^/l) - D^, (55) 



1=0 

where L = 1 + 1/2 and the label 'H' indicates a quantity associated with the metric perturbation in the harmonic 
gauge. To apply this formula in actual calculations (i.e., for a certain orbit on a specific background), one needs to be 
provided with (i) the bare force modes ^pj^^'^^'^'^^^ and (ii) the values of the four regularization parameters Aq, Ba, 
Ca, and Da- To obtain the bare modes ^pj^^^'^^^'^^^ one first calculates the multipole modes of the metric perturbation 
in the harmonic gauge, and then uses the relation 

where k°'l^^^ is the tensor given in Eq. (||) and [/ijj^]^ is the I mode of the trace-reversed metric perturbation (in the 
harmonic gauge) . Whereas this part of the analysis — the derivation of the bare modes — typically involves a numerical 
analysis, the derivation of the regularization parameters may be carried out analytically, bymeans of local analysis 
of the Green's function. This kind of local analysis was described and demonstrated in Sec. 0, where we constructed 
all necessary regularization parameters for two simple cases. 

The first case examined in Sec. ^ was the trivial test case of a static particle in flat spacetime. Here, the mode-sum 
scheme easily reproduced the obvious result of a vanishing self force. We then calculated the regularization parameters 
for the case of a mass particle at a turning point of a radial geodesic in Schwarzschild spacetime. The values of these 
parameters were given in Eq. ( |5^ ) . These values find support from a recent numerical analysis by Lousto |^ 29 1 . 

Our calculation scheme — like the formal prescription by MSTQW on which it relies — is formulated within the 
harmonic gauge. In an accompanying paper pO| | we explore the gauge dependence of the gravitational self force, and 
re-formulate our scheme in terms of other gauges — ones more commonly adopted in perturbation analysis (such as 
the RW gauge). We then conclude that an expression of the form (|5^) is applicable for calculating the self force in 
any specific gauge 'G' (as long as the the self force attains a definite finite value in that gauge), by just replacing 
the harmonic gauge modes ^pj^^^'^'^'^Y^ on the right-hand side of Eq. ( |55| ) with the 'G'-gauge modes [Fa^'^"^'']'^ — the 
ones derived from the G-gauge metric perturbation using Eq. (p6|), with 'H'— > 'G'. The analysis of Ref. [|o) further 
tells us that the regularization parameters in Eqs. ( |55| ) should not carry any gauge label: these parameters are 
"gauge invariant", in the sense that they are always to be derived from the Green's function associated with the 
harmonic-gauge wave operator [the one defined through Eq. (^], irrespective of the gauge in which the bare modes 
are calculated. In particular, we find that the values specified in Eq. (g4|) are valid under any gauge. 



Further applications of the proposed calculation scheme 



The analysis of Sec. Ill provides formal tools for calculating the gravitational self force in any orbit on a Schwarzschild 
background. For any such orbit, the regularization parameters may be derived by local analysis of the Green's function 
modes, based on the separated system of equations ( p^ [supplemented by Eqs. (26)] — as demonstrated in Sec. 0. 



It should be noted, however, that in more general cases than the simple ones considered here, the derivation of the 
regularization parameters shall require much more technical effort: Here, we only had to deal with one equation for 
one Green's function component (in the case of a static particle in flat spacetime), or with three coupled equations 
for three components (in the case of a particle at a turning point of a radial geodesic) . We shall have to face 3 sets of 
7 coupled equations each for a total of 21 components already for an arbitrary point of a radial geodesic. In general 
cases, one would have to deal with up to 58 equations (7 sets of 7 coupled equations each for even perturbation modes, 
and 3 sets of 3 equations each for odd perturbation modes) . 

To fully implement the regularization scheme and derive the self force, one also needs to calculate the bare modes 
of the force. This was already done by Lousto [|l^, using Moncrief's formalism for radial geodesic motion in 
Schwarzschild spacetime. In his analysis, Lousto calculated numerically the (one-sided averaged) values of Aa, Ba, 
and Ca, and presumed a null value for Da based on a proposed zeta- function regularization scheme. To provide a full 
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theoretical basis for Lousto's results in the radial motion case (and verify its consistency with MSTQW regularization) 
we intend to calculate analytically, using our mode-sum scheme, all four parameters Aa, Ba, Ca, and Da for an 
arbitrary point of a radial geodesic Then, the next natural step would be to consider more general, non-radial 
orbits. This would require a considerable amount of both analytic work (deriving the regularization parameters) and 
numerical work (calculating the bare modes of the force by solving the metric perturbation equations). 

We finally comment on the applicability of our mode-sum regularization scheme to orbits in Kerr spacetime. 
Although the theoretical basis for applying our scheme for non-spherically symmetric backgrounds has not been 
developed yet (not even in the scalar case), we believe that a properly generalized version of this scheme could, 
eventually, cope with the Kerr case as well. Such a generalization could still be based on MSTQW's formal prescription 
(^, which applies for any vacuum spacetime. The main obstacle in re-designing our scheme for the Kerr case would 
be, of course, the non-separability of the metric perturbation and Green's function into multipole modes in the time 
domain (such a separation has been a necessary step when executing our scheme in spherically symmetric cases) . This 
difficulty would make both the analytic and numerical parts of the mode-sum scheme more challenging: The bare 
modes of the force would usually have to be calculated in the frequency domain (using the Teukolsky/Chrzanowski 
formalism |28|), and then appropriately summed over Fourier frequencies. As to the analytic part of the scheme, it 
seems to us that, with the use of an appropriate local perturbation analysis, enough information for constructing the 
regularization parameters could be extracted from the time-domain Green's function equations, even without fully 
separating these equations. 
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APPENDIX A: COUPLING OF GREEN'S FUNCTION'S COMPONENTS (SCHWARZSCHILD 

SPACETIME) 

In this appendix we give explicitly the coupling terms appearing in the Green's function equation (|^), for the 
Schwarzschild spacetime case (and using Schwarzschild coordinates). In the following expressions / = 1 — 2M/r and 
a prime denotes d/ dr. 

At^G^, = [r\f'f/'^ - 2/'a,] Gu + 2fdtGtr + [ff" - /(/')V2] G„ + r-^ff {Ggg + sin'^ OG^^) , (Al) 

A'^r'G^, = [/" - r\fT ~ "^r-^f] Gtr + r^f'dtGu + fdtGrr - 2r~^de + cot 9)Gte - 2r-^ sin-^ 9 d^Gt^, (A2) 

^Jf^^G^. = [2f'dr - Af/r^ + /-i(/')V2] Grr + 2f-'f'dtGtr + [f-\f" - /-'(/') V2] Gu 

-4r-^ide + cot 0)Gr9 - 4r-3 sin-^ d^Gr^ + {2r-^ - r-\r' f) {Ggg + sin-^ 9G^^) , (A3) 

A';;G,,, = [f'/r - if + 2f/r)dr - sin-^ e] Gtg + f'dtGre + 2{f /r)deGtr - 2r-2 sin-^ Q cot 9 d^Gt^, (A4) 

KeG^.u = [(/' - 2f/r)dr - 4f/r^ - r'^ sin-^ 9] Grg + 2{f/r)dgGrr + r^'dtGte - 2r-\dg + cot 9)Ggg 

+2r-^ sin"2 q qq^^ _ 2^-2 g q d^Gr^ - 2r-^ sin"^ 9 d^Gg^, (A5) 

■^eeG^.. - h4(//r)5, -I- 2{f/r^) ~ 2r-^ sin-^ e] Ggg + rf-^f'Gtt + {2f - r//')G„ + 4(f/r)deGre 

+2{r-^ cot^ 9 + f'/r) sin^^ 0G^^ - Ar'^ sin^^ 9 cot 9 d^Gg^, (A6) 

Ai^r^G^, - [-Kf/r)dr - 4r-2 cot 9 dg + 2(//r2) + 2r-^ cot^ 9] G^^ + r/- V sin^ 9Gu + {2f - rff) sin^ 9Grr 

+4(//r) sin 9 cos 9Grg + 4(//r)a^G^^ + [2r-^ cos^ 9 + 2{f'/r) sin^ 9] Ggg + Ar'^ cot 9 d^Gg^, (A7) 
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AZG^,, = [^{f + 2f /r)dr ~ 2r-^ coiedg + 2f /r] Gt^ + 2r-^ cot 0d^Gte + 2{f /r)d^Gtr + f'dtGr^, (A8) 

Af^^G^, = [if - 2f/r)dr - 4//r2 - 2r-^ cot 9 de] G,^ + 2(//r)a^G„ + r^f'dtOt^ - 2r-^ sin-^ Q d^G^^ 

+2r-^coted^Gre ~ 2r~^{dg + cot0)Ge^, (A9) 

■^^G^, = [-Kf/r)dr - 2r-2 cot 6 dg ~ 3r-^ cot^ 9 - 2/7r + (2/ - l)/r^] Gg^ + 2{f/r)d^Grg 

+2{J/r){dg ~ 2 cot e)Grv+2r-^ cot ed^Ggg - 2r-^ sibT'^ 9 cot9 d^G^^. (AlO) 



P. A. M. Dirac, Proc. R. Soc. London, Ser. A167, 148 (1938). 

See, e.g., E. E. Flanaga n and R. M. W ald, Phys. Rev. D 54, 6233 (1996); for a pedagogical introduction to the ALD 



equation see E. Poisson, ;r-qc/9912045 



B. S. DeWitt and R. W. Brehme, Ann. Phys. (N.Y.) 9, 220 (1960). 

Laser Interferometer Space Antenna, Proceedings of the Second International LISA Symposium, AIP Conference Proceed- 
ings 456 (American Institute of Physics, Woodbury, 1998). 

See, e.g., E. Poisson and S. W. Leonard, Phys. Rev. D 56, 4789 (1997), and references therein. 
Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D 55, 3457 (1997). 
T. C. Quinn and R. M. Wald, Phys. Rev. D 56, 3381 (1997). 

J. Hadamard, Lectures on Cauchy's Problem in Linear Partial Differential Equations (Dover, New York, 1952). 
T. C. Quinn, Phys. Rev. D 62, 064 029 (2000). 



L. Barack and A. Ori (submitted), gr-qc/0107056 . 

M. J. Pfenning and E. Poisson (submitted), |gr-qc/0012057 

C. O. Lousto, Phys. Rev. Lett. 84, 5251 (2000). 

H. Nakano and M. Sasaki, Prog. Theor. Phys. 105, 197 (2001). 

S. Detweiler, Phys. Rev. Lett. 86, 1931 (2001). 

L. Barack and A. Ori, Phys. Rev. D 61, 061502(R) (2000). 

L. Barack, Phys. Rev. D 62, 084027 (2000). 

L. M. Burko, Class. Quantum Grav. 17, 227 (2000). 

L. M. Burko, Phys. Rev. Lett. 84, 4529 (2000). 

L. Barack and L. M. Burko, Phys. Rev. D62, 084040 (2000). 

L. M. Burko, Y.-T. Liu, and Y. Soen, Phys. Rev. D 63, 024015 (2001). 

L. M. Burko and Y.-T. Liu, Phys. Rev. D 64, 024006 (2001). 

C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San Francisco, 1970). 

Note the mistake in Eqs. (49)-(51) of Ref. where the trace-reversed Green's function has been used instead of the 
Green's function associated with the metric perturbation itself. 

T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). Note that the basis of tensor harmonics introduced by Regge 
and Wheeler is not orthogonal. 

F. J. Zerilh, J. Math. Phys. 11, 2203 (1970); Phys. Rev. D 2, 2141 (1970). 
F. J. Zerilh, Phys. Rev. D 2, 2141 (1970). 

For a comprehensive review of tensor multipole expansions, see K. S. Thorne, Rev. Mod. Phys. 52, 299 (1980). 
P. L. Chrzanowski, Phys. Rev. D 11, 2042 (1975). 

L. Barack and C. O. Lousto, in preparation. Some preliminar y results of thi s work are given in C. O. Lousto, Proceedings 



of the 3rd International LISA Symposium (to be published); gr-qc/0010007 
[30] V. Moncrief, Ann. Phys. (N.Y.) 88, 323 (1974). 



19 



